Neutral atom transport from the termination 

shock to 1 AU 

Maciej Bzowski* and Slawomir Tarnopolski* 

jST , * Space Research Centre PAS, Bartycka 18A, 00-716 Warsaw, Poland 

o 

Cn ' Abstract. Dynamics of H, D, and heavy Energetic Neutral Atoms (ENA) between the termination 

^—^ \ shock and 1 AU is discussed in the context of the forthcoming NASA SMEX mission IBEX. In 

Mh, particular, effects of the velocity-dependent radiation pressure on atomic trajectories are considered 

•^ ' and ionization losses between TS and 1 AU are studied. It is shown, among others, that most of the 

dynamical effects and ionization losses are induced within a few AU from the Sun, which translates 

to the time domain into ^1—3 solar rotations before detection. This loosens considerably time 

requirements for tracking the ionization and radiation pressure history to just prior 3 months. ENA 

seem excellent tracers of the processes within the heliospheric interface, with the transport effects 

^ ■ between the termination shock and detector relatively mild and easy to account for. 

(N 

oo 

1$ INTRODUCTION 

O 

(^ \ This article is a part of reconnaissance for a NASA SMEX mission IBEX (| 1]; McCo- 

■^ ' mas et al, this volume). IBEX will be a spin-stabilized Earth satellite observing Energetic 

P_|. Neutral Atoms (ENA) in the 0.01 - 6 keV energy band, with sensors looking perpen- 
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dicularly to the spin axis directed at the Sun. The prime target of the mission will be 



;-i I imaging of the heliospheric interface via neutral atoms and the goal of an individual 

^ ■ observation will be to obtain the kinematic parameters of the locally registered atom be- 

yond the termination shock (TS). To that end, understanding the processes affecting the 
atoms between their source region and the detector is needed. In particular, a control on 
^ ■ the modifications of the kinematic parameters of the atoms between the source region 

and the detector and a grasp on loss processes underway is necessary. 

Since the ENA gas inside TS is coUisionless at all energies, then - from Boltzmann 
equation - the transport of ENA from TS to the inner heliosphere is governed solely by 
dynamic effects and gain and loss processes. 

The losses include all known ionization processes operating in the heliosphere: charge 
exchange, EUV ionization, and electron impact |2]. For inward-traveling ENA, gains 
(considered by IBEX as unwanted background) include charge exchange between CME, 
CIR and pickup ions and all kinds of neutral populations present inside TS, as well as 
production of ENA in the plasma environment of the spacecraft (McComas et al., this 
volume) and will be neglected. 
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FIGURE 1. Radiation pressure force acting on H and D atoms as function of atomic radial velocity for 
solar minimum and maximum conditions, expressed in the units of solar gravity force ("the jj. units") and 
the net line-integrated solar Lyman-a flux for the past '^ 1 .5 solar cycle, also expressed in the jj. units. 



DYNAMICS OF H AND D ATOMS WITH DOPPLER EFFECT 



From dynamical point of view, heliospheric ENA belong to two categories: sensitive 
and insensitive to solar radiation pressure. The sensitive ones include only hydrogen - 
the most abundant heliospheric element - and its isotope deuterium. All atoms are, of 
course, sensitive to solar gravity. 

All neutral atoms in the heliosphere obey the equation of motion in the form: 



d^r/dt^ = -GM(1 -^{vr,I{t)))r/\r\ 



(1) 



where Vr = {dr/dt) ■r/\r\ is radial velocity and GM the gravity constant times solar 
mass. I(t) is the line-integrated solar Lyman-a flux. For all atoms heavier than D 
the radiation pressure factor jU = and the equation of motion reduces to the simple 
Keplerian case, with the ideal hyperbolic solution. 

The expected range of velocities of heliospheric H and D ENA is ~ 50 — 1000 km/s. 
Beacuse of the geometry of the lines of sight, IBEX is interested only in the inward- 
traveling atoms, with velocities in the source region pointing so that at 1 AU they will be 
at perihelion. Such atoms are being detected at radial velocity with respect to the Sun, 
but cover most of their paths with radial velocities almost equal to their net velocities. 
Consequently, if their net velocity is larger than the spectral range of the solar Lyman-a 
line, the atoms do not experience any radiation pressure at all all the way except within 
a few AU from the Sun, where before detection they reduce their radial velocities to 0. 
Since, however, they are fast, their trajectories will be weakly modified by solar gravity 
and hence close to straight lines. On the other hand, slower ENA (whose velocities 
remain within the spectral range of the solar line), will be affected by radiation pressure, 
which acts against solar gravity. Consequently, they also will follow trajectories close to 
straight lines, except perhaps within a few AU from the detector. 



The behavior of solar Lyman- a line from solar minimum to maximum has been 
recently surveyed by [3]. Based on these data it was possible to fit a general formula 
(Tamopolski, thesis, in preparation): 

IJ,{vrJ) = {a + bI){exp[-cvj] {d + fexp[-gVr-hv^] + jexp[kvr-mv^])) (2) 

which describes the radiation pressure of H or D atoms traveling at radial velocity Vr with 
respect to the Sun. The formula, shown for H and D and solar minimum and maximum 
conditions in the left-hand panel of FigQ is parametrized by the sole parameter /, 
which is the line-integrated solar Lyman-a flux, whose variations during the solar cycle 
(SOLAR 2000, [4]) are shown in the right-hand panel of FigH] 

Due to isotope effect, the deuterium resonant Lyman- a frequency is shifted blueward 
by ~ 80 km/s with respect to H. Hence the D atoms are sensitive to a different part of 
the solar Lyman- a profile than H atoms and due to the mass difference the radiation 
pressure potentially acting on D atoms is additionally reduced two-fold with respect to 
H. The two effects combined will yield different dynamics of the H and D ENA: while 
radiation pressure for H is approximately symmetric with respect to radial velocity, 
this is not the case for D and the effective radiation pressure acting on the incoming D 
atom will be much lower than on a H atom launched from the termination shock with 
the same speed (see Figd left-hand panel). Consequently, the dynamic behavior of the 
incoming D ENA is expected to be more similar to the behavior of purely-Keplerian 
heavy atoms than to that of H atoms. 

Results presented in the following part of the paper are based on numerical solutions 
of Equ.dl]) with the jl factor taken from Equ.©, evaluated for two fixed values of /(?), 
appropriate for solar minimum and solar maximum. For heavy atoms the /i factor was 
set to 0. The atoms were tracked from their perihelions at 1 AU to TS, with velocities 
at perihelion ranging from 50 to 400 km/s. Along with the trajectory tracking the 
probability of ionization for an assumed ionization rate at 1 AU was calculated. 



RESULTS 

Trajectories and energy 

All heliospheric ENA except those from the low end of the IBEX energy band follow 
almost straight line trajectories, with very little velocity change between TS and the 
detector at 1 AU. For detection speeds v >~ 150 km/s there is practically no velocity 
change between TS and 1 AU regardless of the species and solar activity phase (Fig|2l). 
For V lower than this value, H atoms detected at velocities v have TS launch velocities 
vts,h > Vtsd/o, the launch velocities of the heavier species detected at the same v. 
Also, usually v < vts,h, while opposite is true for vxs^d/o- Consequently, mean velocity 
of the H atoms detected at v is higher than the mean velocities of the heavies and their 
travel time is shorter. For example, an O atom detected at 50 km/s needs ~ 22 yr to 
arrive at lAU from the crosswind TS, where it was launched at ~ 26 km/s, while a 
H atom detected at the same 50 km/s needs only ~ 12 yr to make the trip and it was 
launched at ~ 60 km/s. During solar minimum, even the slowest H ENA travel almost 



200 



150 



100 



H, D & O speed at TS vs speed at 1 AU 



00 



70 



50 



30 





























































/ 






























/ 




















































































A 
























































/^ 




























y/ 


/ 


























^ 


/ 


























y^ 


/ 
























y 


y 


/ 
























/" 


/ 






















H, 


max 




/ 


/ 






















H, 


min 


/ 


// 


' 
























D, m 


a^ 


7 
























E 


1, mil 
O 


'/ 

























































Deflection of source position of H, D & O 
at TS from LOS direction 
O 




30 50 70 100 

V at 1 AU (km/s) 



150 200 



25 50 75 100 125 150 175 200 
speed at 1 AU (km/s) 



FIGURE 2. Change of speed of H, D, O atoms at solar minimum and maximum conditions between the 
termination shock and 1 AU (left-hand panel) and angular separation between the line of sight (arrival 
direction at perihelion) and the source position at the termination shock for H, D, O atoms at solar 
minimum and maximum (right-hand panel). The results for O are valid for all heavy atoms. 



constant speed, following straight lines in the whole IBEX energy range. This facilitates 
assessment of the travel times tjs from TS to 1 AU for given detection speed v, which 
can be with good approximation calculated as: 



tjs = tts {B) /v = 125.74/ (1 + 0.253COS 6) /v 



(3) 



where G is offset angle from the nose of TS and for trs in yr v must be in AU/yr. For the 
lowest-energy atoms (~ 100 km/s) the difference in travel time from the nose and tail 
of TS will be about 5 yr, but for typical H ENA at 400 km/s it will be only 1 year. The 
shape of the termination shock was taken from Grzedzielski et al. (this volume). 

Mapping of the launch position of H and D atoms at TS for a given detector pointing 
is practically 1 to 1 , with the deflection angle at the lower end of the energy band within 
~ 10° (Fig|2l right-hand panel). Heavy atoms (insensitive to radiation pressure) are also 
weakly deflected except the slowest ones, which are deflected up to ~ 30° at 50 km/s. 
Deflection of all species drops below 5° for detection speed 100 km/s. This is favorable 
for IBEX since it facilitates linking the pointing of detectors to the source region at TS. 

Also favorable is the relatively short range of the combined solar gravity and solar 
radiation pressure acting on H ENA. The effective "collision time" of H/D ENA with 
the Sun before detection is only ~ 0.25 yr = 3 solar rotations, so the highly fluctuating 
solar Lyman- a output will need to be tracked back no longer than this interval. 



Ionization 



The survival probability of ENA against ionization between TS and 1 AU is a function 
of the launch and detection speeds and of the ionization rate (assumed here to vary with 
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FIGURE 3. Left-hand panel: Survival probability of H (solid lines) and O ENA (broken lines) from 
TS to 1 AU as function of the ionization rate at 1 AU for various detection speeds. The solid lines 
in the left-hand panel correspond to H ENA, the broken lines to heavy atoms (like O); for 200 km/s 
the two lines practically coincide. Right-hand panel: time of effective ionization defined in Equ.(0 for 
H at solar minimum and maximum conditions as function of the ionization rate at 1 AU (right-hand 
panel) for various detection speeds, indicated in the figure, and for solar minimum (solid) and maximum 
conditions(broken). The inset in the right-hand panel shows an example of the course of relative survival 
probability Wrei (t) between TS and r. 



heliocentric distance as 1/r^). For the H ENA expected to be observed by IBEX it is 
quite large in the whole energy band except at the lowest end: for a H atom detected at 
200 km/s even for the highest expected ionization rate 10~^ s~M AU ^ the probability 
is as much as 40%, and for the lowest expected rate equal to 10~^ s^^ it is almost 
1. For faster H, D, and heavy atoms the survival probability is correspondingly larger, 
depending on the ionization rate. On the other hand, the ionization losses will increase 
quickly with the decrease of detection speed, as shown in the left-hand panel of Fig|3l 
In the range of detection speeds shown, one expects some "radiation pressure bonus" 
in the survival probability for the H atoms (and to some extent for D) compared with 
the heavy atoms because of the larger mean speed of H atoms detected at a given v in 
comparison with the mean speed of heavy atoms. This effect is significant only at the 
lowest end of the detection interval. For H, some variations of the survival probability 
with the phase of the solar cycle are expected, but they are not big and probably will be 
masked by fluctuations of the actual ionization rate. 

Most of the ionization losses occur shortly before detection. Denoting by wiau the 
ionization probability of an atom underway from TS to 1 AU and hy w{t) the probability 
of ionization at time t before detection, we can define the relative survival probability 
Wrei (0 ^iid the ionization time Tton as: 



Wrel {t) = l-W (t) /WIAU; W^el {l^ion) =0.8 



(4) 



and determine the ionization times for different species and detection speeds. Results are 
shown in Fig |3] (right-hand panel). The ionization losses suffered by an atom underway 



from TS to 1 AU are relatively small almost all the way down, and increase exponentially 
only within few months before detection. The exact value of t,„„ is a function of 
detection speed and ionization rate, but generally one can say that the 80% ionization 
losses of all atoms occur within 0.1 - 0.25 yr (~ 1 — 3 solar rotations) before detection. 
The ionization times differ most for the lowest expected ionization rates and for the 
lowest energies at detection; for typical H ENA expected from the inner heliosheath the 
ionization time is expected at the level of one solar rotation period. This is favorable for 
IBEX operations because it means that interpretation of a signal requires monitoring of 
the local ionization rate only ~ 1 month backward in time. 



CONCLUSIONS 

The aspects of ENA transport from the termination shock to 1 AU in the contexts of the 
planned IBEX mission can be summarized as follows. 

Both heavy ENA and H and D ENA are excellent tracers of atomic parameters at 
and beyond TS. Velocity change between the termination shock and 1 AU becomes 
noticeable only below 150 km/s and deflection of the trajectories with respect to the 
launch site at TS beyond angular resolution of IBEX instruments for H ENA occurs only 
during solar maximum (which will be beyond IBEX nominal mission time) at the lowest 
limit of the energy band. Heavy atoms are a little stronger deflected, but the deflections 
exceed IBEX resolution only for detection speeds lower than ~ 70 km/s. 

Ionization losses between TS and 1 AU should not be a big problem for H ENA 
expected from the heliosheath. The survival probability drops down below 10% only for 
H ENA traveling at 1 AU slower than 100 km/s. Most of the ionization losses will occur 
within a month or so before detection and hence should be relatively easy to take into 
account. 
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